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We consider a discrete model of population with interaction where 
L ~ — ' the birth and death rates are non linear functions of the population 

,__! size. After proceeding to renormalization of the model parameters, we 

obtain in the limit of large population that the population size evolves 

as a diffusion solution of the SDE 



o 



t ft 



Z? = x+ I f(Z*)ds + 2 I ^ W(ds, du), 
•rf Jo Jo Jo 

o 

en where W(ds, du) is a time space white noise on ([0, oo)) 2 . 

We give a Ray-Knight representation of this diffusion in terms of the 
local times of a reflected Brownian motion H with a drift that depends 

K> upon the local time accumulated by H at its current level, through 

5h the function f'/2. 
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Introduction 

Consider a population evolving in continuous time with m ancestors at time 
t = 0, in which each individual, independently of the others, gives birth 
to children at a constant rate /x, and dies after an exponential time with 
parameter A. For each individual we superimpose additional birth and death 
rates due to interactions with others at a certain rate which depends upon 
the size of the total population. For instance, we might decide that each 
individual dies because of competition at a rate equal to 7 times the number 
of presently alive individuals in the population, which amounts to add a 
global death rate equal to 7(X f m ) 2 , if X™ denotes the total number of alive 
individuals at time t. 

If we consider this population with m = [Nx] ancestors at time t = 0, 
weight each individual with the factor l/N, and choose /i/v = 2iV + 9, X^ = 
2N and 7^ = j/N, then it is shown in Le, Pardoux and Wakolbinger [TT] 
in the above particular case of a quadratic competition term that the "total 
population mass process" converges weakly to the solution of the Feller SDE 
with logistic drift 

dZf = [dZ x t - i{Zf) 2 } dt + 2^/zfdW t , Zl = x. (0.1) 

The diffusion Z x is called Feller diffusion with logistic growth and models the 
evolution of the size of a large population with competition. In this model 9 
represent the supercritical branching parameter while 7 is the rate at which 
each individual is killed by any one of his contemporaneans. This model has 
been studied in Lambert [10], who shows in particular that its extinction 
time is finite almost surely. 

We generalize the logistic model by replacing the quadratic function 9z — 
jz 2 by a more general nonlinear function / of the population size. We then 
obtain in the continuous setting a diffusion which is the solution of the SDE 

ft ft fZ* 

Z? = x+ f(Z*)ds + 2 / W(ds,du), (0.2) 



^0 </o 



where the function / satisfies the following hypothesis. 
Hypothesis A: / e C(R+;R), /(0) = and 3/3 > such that 

f(x + y)- f(x) < /3y Vx,y>0. 



The equation (0.2) has a unique strong solution (see [7]). Note that the 



hypothesis A implies that 

Va;> 0,/(x) < /3x. 
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An equivalent way to write (0.2) is the following. 



Z? = x+ [ f{Z*)ds + 2 f y/Z*dW?, (0.3) 

Jo Jo 

where W x is a standard Brownian motion. However, the joint evolution 
of the various population sizes {Zf,t > 0} corresponding to different initial 
population sizes x would necessitate a complicated description of the joint 



law of the {W x ,x > 0}. Whereas the formulation (0.2) due to Dawson, Li 
[7] with one unique space-time white noise W, describes exactly the joint 
evolution of {Z x ,t > 0, x > 0} which we have in mind. We call this diffusion 
the generalized Feller diffusion. In order to derive this continuous model, 
we first define a discrete model. For defining jointly the discrete model for 
all initial population sizes, we need as in [13] to impose a non symmetric 
competition rule between the individuals, which we will describe in section 
1 below. We do a suitable renormalization of the parameters of the discrete 
model in order to obtain in section 2 a large population limit of our model 
which is a generalized Feller diffusion. Section 3 is devoted to give a Ray 
Knight representation for such a generalized Feller diffusion. The proof of 
this representation uses tools from stochastic analysis, in particular the "ex- 
cursion filtration" , following an analogous proof of another generalized Ray 
Knight theorem in 



1 Discrete model with a general interaction 

In this section we set up a discrete mass continuous time approximation of the 
generalized Feller diffusion. We consider a discrete model of population with 
interaction in which each individual, independently of the others, gives birth 
naturally at rate A, dies naturally at rate /i. Moreover, we suppose that 
each individual gives birth and dies because of interaction with others at 
rates which depend upon the current population size. Moreover, we exclude 
multiple births at any given time and we define the interaction rule through 
a function / which satisfies hypothesis A. 

In order to define our model jointly for all initial sizes, we need to intro- 
duce a non symmetric description of the effect of the interaction as in [3] and 
[TT] . but here we allow the interaction to be favorable to some individuals. 

1.1 The model 

We consider a continuous time Z + -valued population process {X™, t > 0}, 
which starts at time zero from m ancestors who are arranged from left to 
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right, and evolves in continuous time. The left/right order is passed on 
to their offsprings: the daughters are placed on the right of their mothers 
and if at a time t the individual % is located at the left of individual j, 
then all the offsprings of % after time t will be placed on the left of all the 
offsprings of j. Since we have excluded multiple births at any given time, 
this means that the forest of genealogical trees of the population is a planar 
forest of trees, where the ancestor of the population X} is placed on the far 
left, the ancestor of X t 2 — X] immediately on his right, etc... Moreover, we 
draw the genealogical trees in such a way that distinct branches never cross. 
This defines in a non-ambiguous way an order from left to right within the 
population alive at each time t. See Figure 1. We decree that each individual 
feels the interaction with the others placed on his left but not with those on 
his right. Precisely, at any time t, the individual % has an interaction death 
rate equal to (/(Lj(£) + 1) — /(Lj(i))) - or an interaction birth rate equal 
to (/(Lj(i) + 1) — /(Lj(i))) , where Lj(t) denotes the number of individuals 
alive at time t who are located on the left of % in the above planar picture. 
This means that the individual % is under attack by the others located at his 
left if /(Lj(i) + 1) — /(Lj(i)) < while the interaction improve his fertility if 
/(Lj(i) + 1) — f(Li(t)) > 0. Of course, conditionally upon Lj(-), the occurence 
of a "competition death event" or an "interaction birth event" for individual 
% is independent of the other birth/death events and of what happens to the 
other individuals. In order to simplify our formulas, we suppose moreover 
that the first individual in the left/right order has a birth rate equal to 
A + / + (1) and a death rate equal to /i + /~(1). 

The resulting total interaction death and birth rates endured by the pop- 
ulation X t m at time t is then 

VTn vm 

£[(/(*) - f(k - 1)) + - (/(*) - f(k - 1)1 = £(/(*) - f(k - 1)) = f(x?). 

k=l fc=l 

As a result, {X™, t > 0} is a continuous time Z + -valued Markov process, 
which evolves as follows. X™ = m. If X t m = 0, then X™ = for all s > t. 
While at state k > 1, the process 

X m lumps to [ k + 1j at ^ Xk + EJ=i W) - W ~ X )) + ! 
' (k-1, atrate^ + ELiW)-/^" 1 ))"- 

1.2 Coupling over ancestral population size 

The above description specifies the joint evolution of all {X™, t > 0} m >i, 
or in other words of the two-parameter process {X™, t > 0,m > 1}. In 
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L a (t 2 )=9 

Ut,)=4 



Figure 1: Planar forest with five ancestors 
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the case of a linear function /, for each fixed t > 0, {X™, m > 1} is an 
independent increments process. In the case of a nonlinear function /, we 
believe that for t fixed {X™, m > 1} is not a Markov chain. That is to 
say, the conditional law of X r t l+1 given X" differs from its conditional law 
given (X},Xt, . . . ,X™). The intuitive reason for that is that the additional 
information carried by (X^X?, . . . , X t n_1 ) gives us a clue as to the fertility 
or the level of competition that the progeny of the n + 1st ancestor had to 
beneficit or to suffer from, between time and time t. 

However, {X m , m > 1} is a Markov chain with values in the space 
D([0, oo); Z+) of cadlag functions from [0, oo) into Z+, which starts from at 
m — 0. Consequently, in order to describe the law of the whole process, that 
is of the two-parameter process {X™, t > 0, m > 1}, it suffices to describe 
the conditional law of X n , given {X ra_1 }. We now describe that conditional 
law for arbitrary 1 < m < n. Let V t m ' n := X™ — X™, t > 0. Conditionally 
upon {X e , £ < m}, and given that X t m = x(t), t > 0, {V t m ' n , t > 0} is a Z + - 
valued time inhomogeneous Markov process starting from V^ m ' n = n — m, 
whose time-dependent infinitesimal generator {Qk,e(t), k,£ G Z + } is such 
that its off-diagonal terms are given by 

Q 0/ (t) =0, W > 1, and for any k > 1, 

g fc)fc+ i(t) = ^ + ^(/(x(t) + £) - /(x(t) + £ - 1)) + , 
k 

Q k , k -i(t) = Xk + ^2(f(x(t) +£)- f(x(t) + £ - 1))-, 
e=i 
Q M (£)=0, V£^{A;-1,A;,A; + 1}. 

The reader can easily convince himself that this description of the con- 
ditional law of {X™ — X™, t > 0}, given X m is prescribed by what we have 
said above, and that {X m , m > 1} is indeed a Markov chain. 

Remark 1 . 1 iVote that if the function f is increasing on [0, a], a > and 
decreasing on [a, oo), £/ie interaction improves the rate of fertility in a popu- 
lation whose size is smaller than a but for large size the interaction amounts 
to competition within the population. This is reasonable because when the 
population is large, the limitation of resources implies competition within the 
population. For a positive interaction (for moderate population sizes) one can 
realize that an increase in the population size allows a more efficient orga- 
nization of the society, with specalisation among its members, thes resulting 
in better food production, health care, etc... We are mainly interested in the 
model with interaction defined with functions f such that lim^oo f(x) = — oo. 
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Note also that we could have generalized our model to the case /(O) > 0. 
/(O) > would mean an immigration flux. The reader can easily check that 
results in section 2 would still be valid in this case. However in Proposition 



1.3 and in section 3.3 below, assumption /(O) = is crucial, since we need 



the population to get extinct in finite time a.s.. 

1.3 The associated exploration process in the discrete 
model 

The just described reproduction dynamics gives rise to a forest J 7171 of m trees 
of descent, drawn into the plane as sketched in Figure 2. Note also that, with 
the above described construction, the {J rm 1 m > 1), are coupled: the forest 
jrm+i j^g ^g same i aw as thg forest J 7 " 1 to which we add a new tree generated 
by an ancestor placed at the (m + l)st position. If the function / tends to 
— oo and m is large enough, the trees further to the right of the forest J :m 
have a tendency to stay smaller because of the competition : they are "under 
attack" from the trees to their left. From J-" m we read off a continuous and 
piecewise linear M + -valued path H m = (H™) (called the exploration process 
of J 7 " 1 ) which is described as follows. 

Starting from the initial time s = the process H m rises at speed p until 
it hits the top of the first ancestor branch (this is the leaf marked with D 
in Figure 2). There it turns and goes downwards, now at speed —p, until 
arriving at the next branch point (which is B in Figure 2). From there it 
goes upwards into the (yet unexplored) next branch, and proceeds in a similar 
fashion until being back at height 0, which means that the exploration of the 
leftmost tree is completed. Then explore the next tree, and so on. See Figure 
2. 

We define the local time L™(t) accumulated by the process H m at level t 

up to time s by: 

1 f s 
L™(t)=lim-/ l t<H rn <t+e dr. 

The process H m is piecewise linear, continuous with derivative ±p : at any 
time s > 0, the rate of appearance of minima (giving rise to new branches) 
is equal 



PH + 



/(Lf^flr)j + i)-/(Lfwc)j: 



and the rate of appearance of maxima (describing deaths of branches) is 
equal to 

p\ + 



f{[\L^{HT)\+l)~f{[\L^{HT)\) 
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Figure 2: A forest with two trees and its exploration process. 
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Let S m be the time needed in order to explore the forest J 7 " 1 . We have 

S m = inf{s > 0; 1^(0) > m}. 
Under the assumption that S m < oo a.s. for all m > 1, we have the following 



discrete Ray Knight representation (see Figure 1.3) 



(X f m ,£>0,m> 1)= &,(t),f>0,m> 1 



1.4 Renormalized discrete model 



Now we proceed to a renormalization of this model. For x £ K + and N £ N, 
we choose m = \_Nx\ , ft = 2N, A = 2N, we multiply / by N and divide 
by N the argument of the function /. We affect to each individual in the 
population a mass equal to 1/N. Then the total mass process Z N,X , which 
starts from -Lp- at time t = 0, is a Markov process whose evolution can be 
described as follows. 

Z JV ' X jumps from £ to < * ,/ .) /: x , V(- 

J " \^atrate2iVA ; + iVEt 1 ((/(^)-/(^)) > 

Clearly there exist two mutually independent standard Poisson processes Pi 
and P 2 such that 



tN,x 



[Nx\ 1 

nv^ + iv^ 1 




7 iV,2! 



vz,. 
2Ar 2 Z r ^ + N £} (/( 

8=1 



)-/(^r)) 



dr 



N 




NZ' 



2N*Z^ + N £(/(*)- /(^_±))- dr 



j=i 



A' 



Consequently there exists a local martingale M ,x such that 



,JV,i 



N 



+ f f{Z?>*)dr + M?< x . 
Jo 



;i-i) 



Since M '* is a purely discontinous local martingale, its quadratic variation 
[M ,a; ] is given by the sum of the squares of its jumps, i.e. 



[M N > X 



1 

iV 2 " 




NZi 



2N 2 Z^ + N £(/(■£)-/( 



»=i 



"iV" 



f//- 



+A 




VZi. 



2N*Z?* + N ^(/(l)-./(i-Ji))- 



«=i 



iV 



dr 



;i-2) 
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level / 




4>(*) 



i'xpl( iration time s g 1 



12 3 



Figure 3: Discrete Ray Knight representation. 
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We deduce from (1.2) that the conditional quadratic variation (M > x ) of 
M N ' X is given by 

rt ( 1 ^ 

(1.3) 







(M N ' x ) t = 


■J( 


■t 


V 


,x + ^11/11^,0,^'" J 


dr, 


where for 


any z 


_ k 1 _ 

JV> 


V 


Jfc 


ez + 


such that 


k<k', 








ll/lk 


«,«' 


= 


E 


!/(£)-/( 


Z_1 )l- 





i=k+l 



Now we precise the law of the pair (Z N ' X , Z N ' y ), for any x,y E K+ such 
that x < y. Let V ,x ' y := Z ' a — Z N,X , and consider the pair of process 

(Z N > X , V N,x,y \ , which starts from ( -jy-, ^ I at time t = 0, and whose 
dynamic is described by: (Z N ' X , V N ' x,y> ) jumps 



from (A £) to <! 



'(% i) at rate 2iVz + El=i(/(#) " /(¥) + 
;^, i) at rate 2OT + ELi(/(#) - /(¥))" 
- ^ at rate 2iVj + ELi(/(^ ) " /(**=^ + 



,7V' JV / " u J -"" J ^ "-"J ' ^fc=lV V JV J J \ N 

vv]v ,^)atrate2iVj + ELi(/(^)-/( i± F 
The process V N,x ' y can be expressed as follows. 



r N,x,y 



V t N ' x ' y = [Nyl - [iVxJ + f [f(Z?' x + V r N ' x > y ) - f(Z?> x )} dr + M t A 

(1.4) 

where M N,x,y is a local martingale whose conditional quadratic variation 
(M N ' x ' y ) is given by 



1.5) 



Since Z N,X and V N,x,y never jump at the same time, 

[M N ' X , M N ' x ' y ] = 0, hence (M N ' X , M N ' x ' y ) = 0, (1.6) 

which implies that the martingales M N,X and M ,x,y are orthogonal. 
Consequently, Z N > X + V N ' x ' y solves the SDE 



z N,x + yN, X ,y = [Njl + f f ( Z N, X + yN,x,y )dr + j^N,x, y 

JO 



N 
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where M N,x,y is a local martingale with (M N,x,y ) given by 

(M N ' x ' y ) t = (M N ' x ) t + (M N ' x ' y ) t = (M N ' x+y ) t , Vt > 0. 
We then deduce that for any x, y G K+ such x < y, 

z n, x + yl 

In fact, we have that 



7N,x i irN,x,y (fO yN,y 



, /N , x , y _ [Ny\ - [Nx\ 
Vt N 



/ t NV s N ' x ' y , 



£)-/#"+£: 



-±P*\N 



NV, 



N,x.y 



£ ( / W + 4) - /(^ + ' 



fc=i 



TV' 



AT' 



ds 



1 



P 3 I 2iV 2 / V?™ds 



1 



P 4 I 2N 2 / V^ds 



where P 1 , P 2 , P 3 and P 4 are mutually independent standard Poisson pro- 
cesses which are all independent of Z N,X . It follows that conditionally upon 
{Z N,X ,x' < x}, M N,x ' y is a local martingale. 

Remark 1.2 We can also proceed to the renormalisation of the exploration 
process to provide a discrete Ray Knight representation of the process Z N,X . 
We choose the slope p = 2N and we denote by H N the exploration processes 
associated to the forest F N > X of \_Nx\ trees. Let L^(t) be the local time of 
the process H N at level t up to time s. At any time s, the rate of minima of 
H N is equal to 



AN 2 + N 



■ [NL?(H?)\ + _ [NL? (H?)\ 



N 



N 



and the rate of maxima is equal to 



AN 2 + N 



/( mi(HD\ + 1/N) _ f{ lNL?(H?)\ 



N ' ' JK N 

Let S N,X be the time to explore the forest J 7 ^. We have that 

S N > x = M{s>0-L?(0)>^}. 



ds 
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Under the assumption that S N ' X < oo a.s. for all x > 0, the discrete Ray 
Knight representation with the renormalization becomes: 

(z?' x ,t> Q,x > o) = (L% N , x {t),t> 0,x > 0) . 

One could probably deduce from this discrete approximation the Ray Knight 
representation of the general Feller diffusion by a limiting argument, as it is 
done in W in the linear case and in /TT] / in the quadratic case. But in this 
work we use stochastic analysis tools for proving our extended Ray Knight 
theorem. 

1.5 Continous model with a general competition 

Given a space-time white noise W(ds, du), we now define an M + -valued two- 
parameter stochastic process {Z x , t > 0, x > 0} which is such that for each 



fixed x > 0, {Zf,t > 0} is a continuous process, solution of the SDE (0.2) 



We have that for any < x < y, {V t x ' y := Z\ - Zf, t > 0} solves the SDE 

vr = y-x+ \f{z: + v:>y)-f{z:)]ds+2 / w(daM 

Jo Jo Jzf 

(1.7) 

The process V x,y is nonnegative almost surely. We have that f Q f Q s W (ds, du) 
and L f z ° s W(ds, du) are orthogonal since 

[o, z:\ n (z:, z: + v^\ = 0, 

and 

rt ?z% ft f ZS+Vf' v pi r z*+v 3 x ' v 

/ / W(ds,du)+ / W(ds,du)= / W(ds,du) a.s. 

Jo Jo Jo Jz* Jo Jo 

This implies that Z y = Z x + V x ' y a.s. It follows that, for each t > 0, the 
process {Z x ,x > 0} is almost surely non decreasing and for < x < y, the 
conditional law of Z y , given {Zf ,x' < x,t > 0} and Zf = z(t), t > 0, is 



the law of the sum of z plus the solution of (1.7) with Zf replaced by z(t). 
Note that when Z x is replaced by a deterministic trajectory z, the solution 
of (1.7) is independent of {Z x ,x'<x}. Hence the process {Z x , x > 0} 



is a Markov process with values in C([0, oo),R + ), the space of continuous 
functions from [0, oo) into M + , starting from at x = 0. In the case / linear, 
the increments of the mapping x —¥ Z x are independent, for each t > 0. 
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For x > 0, define Tft the extinction time of the process Z x by: 

T x = inf {t > 0; Zf = 0} . 

For any x > 0, we call the process Z x subcritical if it goes extinct almost 
surely in finite time i.e if Tq is finite almost surely. The assumption A implies 
that tW is bounded. Let us introduce the notation 



A(/) := ^"exp (-1 jTM dr j du . (1.8) 

We have the following Proposition. 

Proposition 1.3 Suppose that f satisfies hypothesis A. For any x > 0, Z x 

is subcritical if and only if A(/) = oo. In particular we have: 

i) A sufficient condition for P (Tg < oo) = 1 is: there exists Zq > 1 such 
that f(z) <2,\/z> z , 

ii) A sufficient condition for P (Tq = oo) > is: there exists Zq > 1 and 
5 > such that f(z) > 2 + 5, V z > z . 

Proof: 

Let S G C 2 (1R + ) and < a < x < b. By Ito's formula applied to the 
process Z x and the function S, we have that for any t > 0, 

S(Z x ) = S(x)+ [ {S\Z x s )f{Z x ) + 2S'\Z x s )Z x s )ds + 2 [ S'(Z*)^Z*dW s . 
Jo Jo 

(1.9) 

Let us denote by A the generator of Z x . If we can find a strictly increasing 
function S on the interval [a, b] such that AS = 0, then the drift term in 
(|1.9 ) vanishes and so Z a ' will be just a time changed Brownian motion in 



[S(a), S(b)]. Such a function S is called a scale function of the diffusion Z' 1 
We choose as scale function: for any z > 0, 

S'(z) = / exp \—-z\ dr J du. 

Let us denote by T x the random time at which Z x hits y for the first time. 
We have for any < a < x < b 

HK < I?) = ^ ~ S q [ X \ , and P(r; < oo) = lim P(^' < T x ). 
S(b) — S(a) &->oo 
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If the function S(z) tends to infinity as z goes to infinity, then f{T x < oo) = 
1. Otherwise < P(T^ < oo) < 1. From this we deduce that Z x goes extinct 
almost surely in finite time if and only if lim^oo S(z) = oo i.e. if and only if 
A(/) = oo. The rest of the Proposition is immediate. 



2 Convergence as N — > oo 

The aim of this section is to prove the convergence in law as N — > o o of the 



two-parameter process {Z t ' x , t > 0,x > 0} defined in section 1.4 towards 
the process {Zf, t > 0,x > 0} defined in section 1.5 We need to make 



precise the topology for which this convergence will hold. We note that the 
process Z t ,x (resp. Zf) is a Markov processes indexed by x, with values in 
the space of cadlag (resp. continuous) functions of t D(([0, oo);R + ) (resp. 
C(([0, oo); M+)). So it will be natural to consider a topology of functions of 
x, with values in functions of t. 

For each fixed x, the process t — > Z t ,x is cadlag, constant between its 
jumps, with jumps of size ±iV -1 , while the limit process t — > Zf is continuous. 
On the other hand, both Z t ,x and Zf are discontinuous as functions of 
x. x — > Z x has countably many jumps on any compact interval, but the 
mapping x — > {Zf, t > e}, where e > is arbitrary, has finitely many 
jumps on any compact interval, and it is constant between its jumps. Recall 
that D([0,oo);R + ) equipped with the distance d^ defined by (16.4) in [5] 
is separable and complete, see Theorem 16.3 in |Sj. We have the following 
statement 

Theorem 2.1 Suppose that the Hypothesis A is satisfied. Then as N — > oo ; 

{Z t N ' x , t>0,x>0}^ {Zf, t > 0,x > 0} 

in -D([0, oo); -D([0, oo); M+)), equipped with the Skohorod topology of the space 
of cadlag functions ofx, with values in the Polish space -D([0, oo); M + ) equipped 
with the metric d ^. 

Proof of the theorem 

To prove the theorem, we first show that for fixed x > the sequence 
{Z N ' X , N > 0} is tight in D([0, oo); R+). 
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2.1 Tightness of Z N ' X 

For this end, we first establish a few lemmas. 

Lemma 2.2 For all T > 0, x > 0, there exist a constant Cq > such that 
for all N > 1, 



0<t<T 

Moreover, for all t > 0, N > 1, 



sup E (z t N ' x ) < C . 

<t<T V / 



t 

N,x\ 



E(-y f{Z?*)dr\ < x . 

PROOF: Let (r n ,n > 0) be a sequence of stopping times such that r n tends 
to infinity as n goes to infinity and for any n, I M tA ' x n , t > ) is a martingale 



and ^ iA '^ < W- Taking the expectation on both sides of equation (1.1) at 
time t A r„ , we obtain 



1 n > 



E fe) = 4? 1 + E f P* /(^-)dr ) . (2.1) 



TV " ~\J ■"■'■ 



It follows from the hypothesis Aon/ that 



ei^;„)<^ 



/3 f'mz^jdr 

Jo 



From Gronwall and Fatou Lemmas, we deduce that there exists a constant 
Co > which depends only upon x and T such that 



sup sup E (Z t ' x j < C - 

N>10<t<T \ / 



From (2.1), we deduce that 

-E(f" /( z^),LM. 

Since —f(Z^' x ) > —f3Z^' x , the second statement follows using Fatou's 
Lemma and the first statement. ■ 
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We now have the following lemma. 
Lemma 2.3 For all T > 0, x > 0, there exists a constant C\ > such that 

snpE((M N ' x ) T ) <d. 

N>1 

PROOF: For any N > 1 and k, k! e Z + such that k < k', we set z = ^ and 

ti_ 

N ■ 

k' 



' — We deduce from hypothesis Aon/ that 



N,z,z> — 

i=k+l 



E{(/(^)-/(^i)) +(4)-/<^)) 

= eR/4)-/(^)) + -(/4)-/(^ 1 ) 

i=fc+i L v / \ 

Consequentely, 

|iW < 2£(z / -z)+f(z) -f(z'). (2.2) 



We deduce from (2.2), (1.3) and Lemma 2.2 that 



E ((M^) r ) < £ | (4 + ^)E(Z^) - ^E (/(Z^) } rfr 

<( 4+ ^)C T+-. 
- V N J ° AT 

Hence the lemma. ■ 

It follows from this that M^' 1 is in fact a square integrable martingale. 
We also have 

Lemma 2.4 For all T > 0, x > 0, there exist two constants C2, C3 > such 
that : 



sup sup E ( Z t ,x J < C\ 



V>10<t<T \ Jo 



N>10<t<T 

ft 



sup sup E ( - / Z^ x f{Z^ x )dr ) < C 3 . 



PROOF: We deduce from (1.1) and Ito's formula that 



(Z^ x ) 2 = (JM) + 2 J Z?>*f(Z?>*)dr + (M N >*) t + M?*®, (2.3) 
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where M N ' X '^ is a local martingale. Let (<r n , n > 1) be a sequence of stopping 
times such that lin^oo <r n = +oo and for each n > 1, IM tA ^ , t > 0) is a 
martingale. Taking the expectation on the both sides of (2.3) at time t A a n 



and using hypothesis A, Lemma [2. 3 [ Gronwall and Fatou lemmas we obtain 
that for all T > 0, there exists a constant Ci > such that : 

2 



We also have that 

2E 



sup sup E ( Z t ' x J dr < C 2 . 

7V>10<i<T \ ' 



M " Z^f(Z^)dr) < (Vp- 



C x 



From Hypothesis A, we have —Z^ ,x f{Z^' x ) > —f3(Z^ ,x ) 2 . The result now 
follows from Fatou's Lemma. ■ 

We want to check tightness of the sequence {Z N ' X , N > 0} using Aldous' 
criterion. Let {tn, N > 1} be a sequence of stopping time in [0, T]. We 
deduce from Lemma [2.41 

Proposition 2.5 For any T > and r/, e > ; t/iere exists 5 > suc/i i/iat 

•b+9)AT 



sup sup P 

JV>10<e<<5 



f{Z^ x )dr 



TN 



>V) < e- 



PROOF: Let c be a non negative constant. We have 



(tjv+0)AT 



/(Z^dr 



TJV 



< sup \f(r)\8 

0<r<c 



TiV 



7"JV 



. { ^ >c} |/(Z^)|dr 



But 



"r N +e 



/ hz^ >c} \f{Z?> X )\dr < c- 1 / Z r ^ {f + (Z^ x ) + f-(Z?> x )) dr 

< c- 1 / (2Z^ x f + {Z^ x ) - Z?*f{Z?*j) dr 
Jo 

{2(3{Z^ x f - Z?' x f{Z?> x )) dr. 



<c 



From this and Lemma |2.4l we deduce that V iV > 1 

-{r N +e)/\T 



sup P 

0<6K<5 



TN 



f{Z?' x )dr > 7] < 77-^ 



(rjv+e)AT 



f{Z?")dr 



TN 



< sup UK M + — , 

0<r<c ^ C?7 
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with A = 2/3C2T + C3. The result follows by choosing c = 2A/er], and then 
<J = e77/2sup 0<r<c |/(z)|. ■ 



From Proposition 2.5, the Lebesgue integral term in the right hand side of 



(1.1) satisfies Aldou's condition [A], see pp. The same Proposition, Lemma 



2.2, (1.3) and (2.2) imply that < M N)X > satisfies the same condition, hence 
so does M N ' X , according to Rebolledo's theorem, see [9]. Since all jumps are 
of size jr, tightness follows. We have proved 

Proposition 2.6 For any fixed x > ; the sequence of processes {Z N ' X , N > l} 
is tight in D ( [0, 00) ; K + ) . 



We deduce from Proposition 2.6 the following Corollary. 



Corollary 2.7 For anyO < x < y the sequence of processes {V N ' x ' y , N > l} 
is tight in D ([0, 00); R + ) 

PROOF: For any x fixed the process Z N ' X has jumps equal to ±4 which tends 



to zero as N — > 00. It follows from that and equation (1.1) that any weak 



limit of a converging subsequence of Z N,X is continuous and is the unique weak 
solution of equation ( JQ.2 ) . We deduce that for any x,y > 0, the sequence 
{Z N * - Z N ' X ,N > 1} is tight since {Z N < X ,N> 1} md{Z N > v ,N > 1} are 
tight and both have a continuous limit as iV — > 00. ■ 



2.2 Proof of Theorem 2.1 



From Theorem 13.5 in [5], Theorem 2.1 follows from the two next Proposi- 
tions 



Proposition 2.8 For any n E N, < x\ < x 2 < ■ • • < x r 



(z 



N.xi 



'N,X2 



Z N ' Xn ) ^(Z X1 ,Z 



X2 yx„^ 



as N — >■ 00, for the topology of locally uniform convergence in t. 

Proof: We prove the statement in the case n = 2 only. The general state- 
ment can be proved in a very similar way. For < x\ < x% we consider 
the process ^Z N > X1 ,V N ' X1 ' X2 ), using the notations from section 1 The argu- 
ment preceding the statement of Proposition 2J3 implies that the sequences 
of martingales M N,Xl and M N,Xl,X2 are tight. Hence 



(Z N >*i,V N ' xl > X2 ,M N > xl ,M N > x i> x *) is tight. Thanks to JTT}, JO(, JO) , ( [To] 
and (1.6), any converging subsequence of 



2 CONVERGENCE AS N -> oo 



20 



[Z N >*>- 1 V N >* 1 > xa ,M N < sl ,M N < xl * 2 ,N> 1} has a weak limit 
(Z x \V Xl ' x \M x \M Xl ' X2 ) which satisfies 

Z^=x 1+ f /(Zf)ds + Mf 
Jo 

V™ =x 2 - Xl + f f [f{Zf + V?**) - /(Zf )] ds + Mf 1 ' xa , 
where the continuous martingales M X1 and M Xl,X2 satisfy 

(M x ) t = A f Z Xl ds, (M Xl ' X2 ) t = A f V XuX2 ds, (M x \M xl ' X2 ) t = 0. 
Jo Jo 

This implies that the pair (Z Xl , V X1,X2 ) is a weak solution of the system of 



SDEs (0.2) and (1.7), driven by the same space-time white noise. The result 



follows from the uniqueness of the system, see [7j . ■ 

Proposition 2.9 There exists a constant C , which depends only upon 9 and 
T, such that for any < x < y < z, which are such thaty — x < 1, z — y < 1, 



E 



sup \Z t - Z t | x sup \Z t - Z t 

0<t<T 0<t<T 



< C\z-x\ 



We first prove the 
Lemma 2.10 For any < x < y, we have 



sup E (z? ' y - Z. 

<t<T V 



0<i<T 



< )= o^ W -KM N 



f [Ny\ [Nx\ 



N,x,y^ < | L Jt ^J _ L JV -^J ] e /3T 



PROOF: Let (a n , n > 0) be a sequence of stopping times such that lim^oo a n = 
+oo and ( M tA ^ y ) is a martingale. Taking the expectation on the both sides 
of (1.4) at time t A a n we obtain that 



Using Gronwall and Fatou lemmas, we obtain that 



(2.4) 



sup E(VJ 

0<t<T 



N,x,y\ 



< 



f[Ny\ [Nx\ 



\ N 



N 



JT 
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Proof of Proposition 2.9 Using equation (1.4), a stopping time argu- 



ment as above, Lemma [2 .10| and Fatou's lemma, where we take advantage of 
the inequality f{Z^ x ) - f(Z?' x + V r N > x > y ) > -/3V r N ' x ' y , we deduce that 



E 



jT [f(Z?*) - f{Z^ x + V r N **)] dr^j < 



[Ny\ [Ny\ 



N 



N 



(2.5) 



We now deduce from (1.5), Lemma 2.10, inequalities (2.5) and (2.2) that for 



each t > 0, there exists a constant C(t) > such that 

E «M»«)«)<c W (!^J-!f! 



(2.6) 



This implies that M N,X,V is in fact a square integrable martingale. For any 



N,z 



< x < y < z, we have Z t ' — Z 



?N,y 



v; 



N,y,z 



N,y 



and Z t — Z t 



zAf.a: 



K 



JV,!B,J/ 



for 



any t > 0. On the other hand we deduce from (1.4) and the hypothesis A 

2 



sup (V^) 2 < 3 (^ll - i^i) + 3/3 2 T / sup (Vf™) 2 dr 

<t<T \ N N ) JO 0<s<r 

+ 3 sup (m?* 



sup {V t N n 2 < 3 (^P- - ^f) + 3/3 2 T / sup {V^fdr 

<t<T V iV iV / -/O 0<s<r 



+ 3 sup (M t N ' y ' z 



o<t<r 



Now let Q x,y := a (Z t ,x , Z t ' y , t > J be the filtration generated by Z N ' X and 

Z^. It is clear that for any t, V t ,x ' y is measurable with respect to Q x,y . We 
then have 



E 



sup I V t " I x sup I V t 

0<t<T 0<t<T 



sup |Vj y | El sup |v/ t ' w ' \\y ,y 

0<t<T \0<t<T 



Conditionally upon Z N,x and Z N,y = u(.), V N ' y,z solves the following SDE 



v; 



N,y,z 



[Nz\ - [Ny\ 

N 



+ / [f(V r N ' y > z + u(r))-f(u(r))}dr + M t N > y > z , 



where M N,y,z is a martingale conditionally upon Q x,y , hence the arguments 
used in Lemma [2.101 lead to 



l(l ^,),(lM_l&l) 



3 /3T 
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and those used to prove (2.5) yield 



E 



/ f(Z?> v ) - f(Z?' y + V r N ' y > z )dr\g x A < 



[Nz\ [Ny\ 



N 



N 



From this we deduce (see the proof of (2.6)) that 



e ((M N >y>*) t \g*>y) < c(t) ( iM - iM) 



N 



From Doobs's inequality we have 



e( sup \M t N ' y ' z \ 2 \g x A =E((M N ' y ' z ) T \g x ' y ) 

\0<t<T J 



Since < z — y <1, we deduce that 



e (^ i^Tie-) < 3d + c(D) (i^i - IM) 

+ 3/3 2 T f E ( sup (y/to*) 2 ^** ) dr, 

JO \0<s<r / 

From this and Gronwall's lemma we deduce that there exists a constant 
K\ > such that 



E I sup \V t N ' y ' z \ 2 \g x ' y 1 < 



0<t<T 



f[Nz\ [Ny\\ 
Kl \~N N~)- 



(2.7) 



sup (V?™) 2 

0<t<T 



<K ( [Ny\ [Nx\ 
~ Kl I ~N N~ 



Similary we have 
E 
Since < y — x < z — x and < z — y < z — x, we deduce that 

E 



sup |Vj S | x sup VJ y | 



_0<t<T 

hence the result. 



0<t<T 



< ^ 2 / IA*J L^J 



iv 



A' 



Proof of Theorem 2.1 We now show that for any T > 0, 



tN,X 



{Zp x , < i < T, x > 0} =► {£*, < £ < T, x > 0} 
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in £)([0, oo); .D([0, T],M+)). From Theorems 13.1 and 16.8 in [5J, since from 



Proposition |2.8[ for all n > 1, < X\ < ■ ■ ■ < x n , 

{Z N ' X \. . . , Z N ' Xn ) =► {Z X1 , ...,Z Xn ) 

in D([Q, T];R n ), it suffices to show that for all x > 0, e, r] > 0, there exists 
JV > 1 and S > such that for all JV > JV , 

%,j(^) > e) < 77, (2-8) 

where for a function (x, t) — > ^(x, t) 

%,,5(2) = sup inf {||2;(x, •) - z(xi, -)||, lk(^2, •) - ^, Oil}; 

with the notation ||^(x, -)|| = sup 0<t<T |z(x,t)|. But from the proof of Theo- 
rem 13.5 in [3], (2.8) for Z N follows from Proposition 2.9 



3 Ray Knight representation of a general Feller 
diffusion 

In this section we establish a Ray-Knight representation of Feller's branching 



diffusion solution of (0.2), in terms of the local time of a reflected Brownian 
motion H with a drift that depends upon the local time accumulated by H 
at its current level, through the function /' where / is a function satisfying 
the following hypothesis. 

Hypothesis B: / G C 1 (M+; K+), /(0) = and there exist a constant (3 > 
such that 

fix) <fi Vx > 0. 

Note that hypothesis B follows from hypothesis A if we assume that / is 
differentiate. 

The proof we give here is purely in terms of stochastic analysis, and is 
inspired by previous work of Norris, Rogers and Williams [12] and Pardoux, 
Wakolbinger [T3]. We specify an SDE for a process (H s ), from which the 



generalized Feller's diffusion solution of (0.2) can be read off from reflected 
Brownian motion with a drift that is a function of the local time accumulated 
by H at its current level. One way to understand the form of the drift is to 
see (H s ) as the limit of the exploration process H N of the forest of random 
trees associated to Z N)X . Precisely, fix z G C(R + ;R + ), the set of continuous 
functions from M. + into M. + and consider the stochastic differential equation 

HI = B * + \[ f( z ( H r) + L r( H r))dr + ^(0), (3.1) 
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where B is a standard Brownian motion, and for s, t > L z s (t) is the local 
time accumulated by H z at level t up to time s. For x > define 

S^ = inf {r > : L z s (0) > x} and S = sup S" x . 

We first suppose that / satisfies hypothesis B and the following. 
Hypothesis C: 

3 a, beR:Vz>0, \f'{z)\ < az + 6. 

3.1 Case where /' satisfies hypothesis C. 

In this subsection we suppose that / verifies hypothesis C. We have 



Proposition 3.1 For any z e C(R + ;R + ), equation (3.1) has a unique weak 
solution. 

PROOF: Let H denote Brownian motion reflected above defined on a 
probability space (Q, J 7 , P). H solves the following equation 

H S = B S + h s (0), 

where 5 is a J s standard Brownian motion, and L is the local time of H. 

Let 

M s-=\j f'(z(H r ) + L r (H r ))dB r and G s = exp (m s - \{M)}j . 

We will show below that E(G S ) = 1, for all s > 0, which is a sufficient 
condition for G to be a martingale. By application of the Girsanov theorem, 
there exists a new probability P 2 on (Q, J 7 ) such that 

df\ 

^l* = G..«>0, 

where (J-" s , s > 0) is the natural filtration of H. Moreover under P 2 , 

B\ := B s - \ [ f'(z(H r ) + L r (H r ))dr, s > 
^ Jo 

is a standard Brownian motion. The fact that E(G S ) = 1 for any s > 
follows thanks to assumption C from the existence of a constant c such that 

sup E(exp(c(L r (iJ r .)) 2 ) < oo. (3.2) 

0<r<s 



The inequality (3.2) is estabilished in [13], see Lemma 2 and Lemma 3. The 
uniqueness is also proved in [T3] and that argument does not make use of 
hypothesis C. ■ 
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For K > 0, we now consider Brownian motion reflected in the interval 

[0, K] 

H K = B a + ±Lf{0)-±L?(K-), 

defined on (O, J 7 , P), where L K denotes the local time of H . Define for 
x > 

Sf = inf (s>0;Lf(0) > x) . 

We again define 

Mf := i jf /'(*(#*) + L?(H r K ))dB r and Gf = exp (m? - \{M K ))j . 

The same argument as above shows that K(G^) = 1 for all s > 0. This 
implies that there exists a probability F K,Z defined on the measurable space 
(ft, J 7 ) under which 

Bl = H K - hf (0) + iif (IT) - i J' f'{z{H?) + L?(H r K ))dr, 

is a P^-Brownian motion. That is the equation 

H K = B S + hf (0) - iLf (K-) + ± jf' /'(*(#*) + L«(H«))dr (3.3) 



admits a weak solution. Uniqueness of the weak solution of (3.3) is obtained 



in a similar way as concerning (3.1). Moreover we have that (see again |13j ) 

F K > Z (S? < oo) = 1. 
We now have the following Ray Knight representation. 

Proposition 3.2 For any K > and z e C(R + ; R + ), the law of < L^ K (t), 

< t < K > under F K ' Z is the same as the law of {Z^' z , < t < K}, where 
Z x ' z solves the SDE 

dZt' z = [f(Zr + z{t)) - f(z(t))} dt + 2y/%*dW t , Zr = x, (3.4) 

and W is a standard Brownian motion. 

The proof of this Proposition is similar to that done in [13] in the quadratic 
case. We will give below some details of the proof in the more general case 
without hypothesis C. 
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3.2 Existence and uniqueness of weak solution of (3.3) 
without hypothesis C 

Now we do not suppose anymore that / satisfies hypothesis C. We have 

Proposition 3.3 For any K > 0, z G C(M + ; M. + ), there exists a probability 



ys>K,z un d er w fii c fi equation (3.3) has a unique weak solution on the random 
interval [0,S K ), where S K = sup x>0 S^. 

PROOF: Consider again, for K > 0, the Brownian motion H reflected in 
the interval [0,K] defined on (fi, J*,P). 

H K = B s + 1 -L?(0)- 1 -L?(K-). 

For n > 1, we define the function g n (r) = f'(nAr). It is clear that there exist 



two constants a,b > such that 1^(^)1 < ar + b. Thanks to Proposition 3.1 
there exits for each n > 1 a probability p-^- 2 '™ such, 

jjn>K,z,n ( i ^ 

<Al k=exp \M*»-hM*») t \,s>0, 



dP ' s ^ t 2 

where M^ n := \ £ g n (z(H?) + L?{H*))dB r , Under P^- n , 

4^ = #f - \lf (0) - ^ j\n{z{H«) + L?(H r K ))dr + h?(K-),Vs > 0, 

(3.5) 
is a standard Brownian motian. For n > 1, we define the stopping time 

T n = inf i s > 0; sup [z(t) + Lf (£)] > n I . 

We need the following result which is a variant of Theorem 1.3.5 from Stroock- 
Varadhan [15], whose proof is very similar to that in [15J. 

Theorem 3.4 Let Q = C(R+,R+) be the canonical path space with its 
canonical filtration {Ft}, and let (T n ) be an increasing sequence of stopping 
times satisfying T n < S K a.s. Vn > 1. Suppose there is a sequence (P n ) of 
probabilities on (Q, J 7 ) such that 

• P„+i agrees with P n on J-r n ; 

• for each x > 0, 

P n (T n < S? ) -> as n -»- oo. 
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Then there exists a probability P on (fl, J~sk) such that for each n, 

P = P n on 7 Tn . 

This proves the existence of a probability F K ' Z on (Q, J-gx), provided we show 
that for all x > 0, 

F K ' z ' n (T n < Sf ) -> as n -> oo. 

We have 

^' 2 ' n (T n < Sf ) = P x ' 2 ' n ( sup L* K (t) > n ) . 



P' 



From Propostion 



3.2 



;2^' B ) < t < /(Tsolution of the SDE 



under f> K > z ' n ^ (L^ K (t), < t < K) has the same law as 



/ rz(t>+z, t \ 

dZ x t ^ n =\ g n {u)du)dt + 2^/zf^dW t , Zl 

\Jz(t) I 



x,n 

- X. 



For any x > 0, consider the process Z x , which is solution of the SDE 



ft rt 

Z? = x + l3j Z*dr + 2 \lZ;:J\\'r. 



By a well known comparison theorem for one dimensional SDEs, see 
theorem X.3.7, we have that for any x > and z G C(K + ; K+), Z x,z,n < Z 
a.s. . We then have that 

$K,z,n ( Tn < S K\ = ^K,z,n f gup jK^ > \ 

\0<t<K x J 



P ( sup Z t ' ' > n 

yO<t<K 



< P sup Zf > n\ -)• as n -)• oo. 

\0<t<K ' J 

We thus have proved for all K > and z G C(R + ;R + ), the existence of a 
probability F K,Z under which on [0, S ), 

SJ = #f - \lf (0) + ± jT /'(^ r A ') + L«{H?))dr + iif (K-) 

is a standard Brownian motion. Uniqueness is obtained in a similar way as 
in [13]. Hence the Proposition. ■ 
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For any z G C(R + ; K+), we have the following Ray Knight representation. 

Proposition 3.5 For any K > 0, z G C(R + ;R + ) and x > 0, the law of 
(^fir(t), < t < K\ under F K > Z is the same as the law of {Z^\ < t < K) 

PROOF: For K > and z G C(K+), we work under the probability measure 
F K ' Z . Using Tanaka's formula, we have for any r > and < t < K, the 
following identity 

(H? - t)- = (-t)~ + jf" ! {H K< t} dH? + l -L? (t) (3.6) 



Recall that for any x > 0, F K > Z (S£ < oo) = 1. Hence from gg , 



L*K{t) = 2 / ! {H K< t} dHf. 

Jo 



Combining with equation (3.3), we get 



fS K rS K 



L^(t)=x + 2 [ * l {H K <t} dB s + [ * l {H K <t} f{z{Hf)+L«{H?))ds. 
Jo Jo 

From the generalized occupation time formula (see Exercise ^1.1.15 in 
we obtain 



'"O™ /* C /*o. 



/ ' l {HF< Qf'{z{H? ) + Lf (tff ))<fe = / / " /' (*(r) + Lf {r))dLf (r)dr 
Jo Jo Jo 



t 

K 



f(z(r) + L««(r))-f(z(r)))dr. 

o 

The key idea of the proof is now to introduce the "excursion filtration" , as 
in [12] and [13]. For any < t < K and s > 0, let define 



A s (t) := / l {H K<t}dr, r{r,t) = inf {s > 0;A s {t) > r} 
Jo 

J{s,t) := f l {H? < t} dB r £(r,t) := J(r(r,t),t) 
Jo 

F(r,t) '■= <r{£{r,t) ■ < u < r), q t = JF^t), 



N t '■= / l{H«<t}dB s . 
Jo 

For fixed t, the process (J(s,t),s > 0) is a martingale with respect to J 7 ,, 
while (£(r, t),r > 0) is a J-^^-martingale and its quadratric variation equals 
r. Consequently (£(r, t),r > 0) is a J-"(r, t)-Brownian motion. We then have 
the 
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Lemma 3.6 The process (N t , <t < K) is a continuous St-martingale with 
its quadratic variation given by 



(N)t = A [ L« K (r)dr. 
Jo 



This result is an easy consequence of Theorem 1 in [12] . By the martingale 
representation theorem, we deduce that there exists a Brownian motion W 
such that 



'o 
Consequently for all < t < K, 



N t = 2 f ^L^(r)dW r . 



L%(t)=x + j (f{z{r) + L K s x{r))-f{z{r)))dr + 2j ^L««(r)dW r . 



3.3 Ray Knight theorem in the subcritical case 



We first prove the following proposition (recall the definition (1.8) of A(/)) 



Proposition 3.7 Suppose that f satisfies hypothesis B and A(/) = oo. 



Then equation (3.1) admits a unique weak solution on [0,5*). 



PROOF: For x > and K > 0, let define 

n K >* = J sup H K ' <K, V K' > K 

{ [o,&f ] 

For any x > and z G C, since we are in the subcritical case, there exists 
T XjZ < oo a.s. such that , Z*' z — 0, V t > T XjZ . We deduce from this and 



Proposition |3.5| that for any fixed x > 0, 

Q = U K>0 Q K ' X a.s. . 

Note that the family of events {tt K ' x , K > 0} is increasing, and on Q K ' X , 
H = H a.s., for any K' > K. We can define a probability F z ' x on (0,,J r s x ) 
such that p^^=p^ 2 on Vt K > x . Under P 2 ' x , on [0, S x ] 

B\ = H S -IJ S f'(z(r) + L r (H r ))dr - ^L s (0) 
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is a standard Brownian motion. This proves that (3.1) has a weak solution 
on [0, S x ] whose uniqueness can be proved as in [13J. ON A A NOUVEAU 
BESOIN DU TH DE SV ! We can deduce that there exists a probabilty P z 
under which, on [0, S) B z is a standard Brownian motion, where 

S = supS^. 

x>0 



For z = 0, we write P = P 2 . The following statement is a generalized Ray 
Knight theorem in the subcritical case. 

Theorem 3.8 Suppose that f satisfies Hypothesis B and A(/) = oo. Then 
the law of the random fields {L Sx (t),t > 0,x > 0} under the probability P is 
the same as the law of{Zf,t >0,x> 0}. 

We first establish the following Proposition. 



Proposition 3.9 Assume that the two assumptions of Theorem \3.8 holds. 



Then for any x and z G C(R + ; K + ) fixed, the law of {L Sx (t),t > 0} under P 2 
coincides with is the law of {Z^ ,z ,t > 0}. 

Proof: We have that for any K > 0, z E C(R + ; R+), under P^' 2 , (L^ K (t), 

< t < Kj has the same law as (Z^ ,z , < t < K). A consequence of this is 
that for any < K < K', 

{L$K(t),0 < t < K} ( = } {Lf#{t),0 < t < KJ. (3.7) 

It now follows that for any K, under P 2 , (Ls x (t),0 <t< K) has the same 
law as (Lg K (t), < t < K J under P x ' 2 ' x . We then obtain that for any K > 

(L Sx (t),0<t<K) { ^(Z^,0<t<K). 

Hence the proposition, letting A' go to oo. ■ 

In particular, for x fixed, the law of {Ls x {t),t > 0} under P is the same 
as the law of {Zf, t > 0}. 



Remark 3.10 The identity (3.7) could also obtained from a generalization 
of Lemma 2. 1 in Delmas $j$. For < a < b, we define the application Tc a ' b 
with maps continuous trajectories with value in [0,6] into trajectories with 
values in [0,a] as follows. If u G C(R + , [0,6]), 

p u {s) = l u (r)<adr, 7t a ' b (u)(s) =u(p~ 1 (s)). 
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The following equality in law holds. 

ir a >\H b ) ( = } H a . 
This identity together with the strong Markov property of the Brownian mo- 



tion implies (3.7). 



Proof of Theorem 3.8 Recall that (Z x ,x > 0) is a Markov process with 



value in the space of continuous paths from R + into R+ with compact sup- 



port. From Proposition 3.9 with z = 0, its marginal laws coincide with 
those of Ls x (-)- We now check that (Lg x (.), x > 0) is a Markov process. This 
follows readily from the fact that for any < x < y, conditionnaly upon 
(Ls x ,(.),x' < x) and given L Sx (.) = z(.), on [0, S y ] the process H x s := H Sx+s 
solves the SDE 

H x s =Bs + \f Q U'{z{H x r ) + L" r {H?))) dr + llftO), 

where B is a Brownian motion independent of (Ls ,{t),x' ^ x,0 <t < S x ) 
and L z denotes the local time of H x , which is also the additional local time 
accumulated by H after time S x . To complete the proof of the theorem it now 
suffices to prove that for any x, y > the conditional law of [L Sx+y {t),t > 0) 
given (Ls x (t),t > 0) is the same as the conditional law of (Z x+V ,t >) given 
(Z x ,t > 0). Conditioned upon Ls x (.) = z(.), Ls x+y (-)—Ls x () is the collection 
of local times accumulated by H x up to time S y , and it has the same law as 
L s (.) while conditionally upon Z x = z(.), Z x+y — Z x has the same law as 



Z y,z . The identity of those two laws has been established in Proposition 3.9 



We can deduce from the Proposition 3^ and the occupation time formula 

that 

Corollary 3.11 Suppose that f satisfies Hypothesis B and A(/) = oo. We 
have 

Va;>0, f{S x <oo) = l. 

PROOF: Let g(h) = 1, for any h > 0. By the occupation times's formula, we 
have 

S x = g(H r )dr 

Jo 

"OO 



/ L Sx {t)dt = / Z x dr < oo a.s. 
Jo Jo 
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Note that S x is the total mass of the process (Zf,t > 0). 

References 

[1] Aldous, D. (1978). Stopping times and tightness. Ann. Probab. 6, 335- 
340. 

[2] Aldous, D. (1991). The continuum random tree I. Ann. Probab. 19, 
1-28. 

[3] BaM., Pardoux E. (2012). The effect of competition on the height 
and length of the forest of genealogical trees of a large population, to 
appear. 

[4] Ba, M., Pardoux, E. and Sow, A.B. (2012). Binary trees, explo- 
ration processes, and an extended Ray-Knight Theorem. J. App. Probab. 
49, 201-216. 

[5] Billingsley, P. (1999). Convergence of Probability Measures, 2nd ed. 
John Wiley, New York. 

[6] Cattiaux P., Collet P., Lambert A., Martinez S., Meleard 
S., San Martin J. (2009) Quasi-stationary distributions and diffusion 
models in population dynamics, Ann. Probab. 37, 1926-1969. 

[7] Dawson D.A, Li, Z.H. (2012) Stochastic equations, flows and 
measure- valued processes., Ann. Probab. 40, 813-857. 

[8] Delmas, J.F. (2008) Height process for super-critical continuous state 
branching process, Markov Process. Related Fields 14, 309-326. 

[9] Joffe, A., Metivier, M. (1986). Weak convergence of sequences of 
semi-martingales with applications to multitype branching processes. 
Ad. Appl. Probab. 18, 20-65. 

[10] Lambert, A. (2005). The branching process with logistic growth. Ann. 
Probab. 15, 1506-1535. 

[11] Le, V., Pardoux, E. and Wakolbinger, A. (2012). Trees under 
attack: a Ray-Knight representation of Feller's branching diffusion with 
logistic growth. Probab. Theory & Relat. Fields, to appear. 



REFERENCES 33 

[12] Norris, J.R., Rogers, L.C.G. and Williams, D. (1987). Self- 
avoiding random walks: a Brownian motion model with local time 
drift. Probab. Theory & Relat. Fields 74, 271-287. 

[13] Pardoux, E. and Wakolbinger, A. (2011). From Brownian motion 
with a local time drift to Feller's branching diffusion with logistic growth. 
Elec. Comm. in Probab. 16, 720-731. 

[14] Revuz, D., Yor, M. (1999) Continuous Martingales and Brownien 
Motion, 3d ed., Grundlehren der mathematischen Wissenschaften 293, 
Springer Verlag. 

[15] Stroock, D.W. and Varadhan, S.R.S. (1979). Multidimensional 
diffusion processes, Grundlehren der mathematischen Wissenschaften 
233, Springer- Verlag. 



